Intermittency in a single event 
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Abstract 

The possibility to study intermittency in a single event of high multiplicity 
is investigated in the framework of the a— model. It is found that, for cascade 
long enough, the dispersion of intermittency exponents obtained from individual 
events is fairly small. This fact opens the possibility to study the distribution 
of the intermittency parameters characterizing the cascades seen (by observing 
intermittency) in particle spectra. 
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1. The original suggestion of intermittent behaviour in multiparticle pro- 
duction at high energies [|IJ was based on analysis of a single event of very 
high multiplicity recorded by the JACEE collaboration ||. It was soon re- 
alized, however, that the idea can be applied to events of any multiplicity 
provided that a proper averaging of the distributions is performed ||. This 
led to many successful experimental studies of intermittency , and allowed 
to express the effect in terms of the multiparticle correlation functions |J. It 
should be realized, however, that the averaging procedure, apart from clear 
advantages, brings also a danger of overlooking some interesting effects if 
they are present only in a part of events produced in high-energy collisions. 
It seems therefore interesting to study intermittency parameters of individual 
events, hoping that they may indicate some specific production mechanism (a 
typical example is the production of quark-gluon plasma which is expected 
to be characterized by specific intermittency exponents, see e.g. 0, and 
certainly not expected to be present in each event). 

Such studies should necessarily be restricted to high-multiplicity events 
because only there one may expect the statistical fluctuations to be under 
control. However, even neglecting statistical errors due to the finite num- 
ber of particles, there remains an intrinsic uncertainty of the intermittency 
parameters: the cascade responsible for intermitent behaviour has different 
realizations in different events. As the intermittency exponents determined 
from different realizations of the same random cascade are expected to scat- 
ter around the average, the method has a finite resolution with respect to 
the parameters of the random cascade. Clearly, the resolution is a function 
of the number of steps in the cascade. 

In the present paper we investigate the distribution of intermittency ex- 
ponents obtained from analysis of individual events, using as a tool the 
a— model of one- dimensional random cascade Pj. We concentrate on two 
problems : 

(a) how much the average value of an intermittency exponent obtained 
from analysis of individual events differs from its "theoretical" value cal- 
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culated from the assumed parameters of the random cascade and from the 
"standard" value obtained by averaging factorial moments over many events. 

(b) what is the dispersion of this distribution or, in other words, what 
is the resolution of the measurement and how it depends on the number of 
steps in the cascade. 

We have investigated cascades of up to 10 steps. Our results are described 
in some detail below. They can be summarized as follows : 

(i) the average value of the intermittency exponent reproduces well the 
value obtained by averaging the factorial moments over events. This result 
weakly depends on the number of steps in the cascade. However, both num- 
bers have a tendency to underestimate the theoretical value for cascades of 
short lengths. 

(ii) the dispersion of the distribution is inversely proportional to the 
length of the cascade. For cascades of less than 6 steps it becomes fairly 
large, but still substantially smaller than the average value. 

We thus find that a measurement of the distribution of intermittency ex- 
ponents obtained from individual events of high multiplicity may indeed help 
to reveal existence of groups of events emerging from cascades with different 
characteristics. 

2. The a— model of random cascading |l| describes a multiparticle event 
as a series of steps in which each phase-space interval is divided into some 
number of equal parts. At any step n ( n — 1, 2, . . . , N) particle density in 
each of the parts is obtained by multiplication of the density at the (n-l)th 
step by one of the two values (a,b) of random variable W with the probabil- 
ities a and /3, respectively. For simplicity one assumes also : 



where <> denotes the average value henceforth. Note that (1) implies : 



< W >= aa + (3b = 1 
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So that the model is defined by two parameters a and b. 

In our simulation of the a-model we have divided each bin into 2 parts, 
so the number of bins at each step equals : 

M(n) = 2 n (3) 

and thus the length of each bin is equal to : 

where D is the total phase-space interval. 

The "standard" method is to study scaling behaviour of the normalized 
moments of particle densities : 

< Zl{d) >=< (x m (d)) q > (5) 

Here x m (d) is the density obtained after n steps of the cascade in the mth 
bin ( m = 1, . . . , M(n) ), and the average is taken over all considered events. 
It follows from (1) that < >= 1. In the a— model < Z^id) > follows 
the power law : 

2D 

< zm >= (^-r (6) 

where the intermittency exponents are given by : 

ip q = log 2 < W q > (7) 



3. If one is interested in event-by-event analysis, one is forced to consider 
the so-called horizontal average Z Q (d) [|TJ : 

M 

Z\d) = M" 1 £ Zl(d) (8) 

m=l 

obtained by averaging over all bins. In the a— model the average particle 
density is independent of m and thus Z q (d) follows the same scaling law (6) 

as Z^(d) : 

2D 

Z%d) = (— r (9) 
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As we have already explained, ip q calculated from (9) fluctuate from event 
to event even for fixed parameters of the cascade. Its average over many 
events should approach the value given by (7). The dispersion around the 
average, however, does not vanish, even in the limit of infinite number of 
events. In other words, even for events with very large multiplicity we cannot 
determine intermittency exponents with arbitrarily high precision: there is a 
"natural" uncertainty of this measurement. This uncertainty is expected to 
decrease with increasing number of steps in the cascade. Furthermore, the 
dispersion of the distribution of the factorial moment Z q (d) can be estimated 
as: 

D 2 (Z q (d)) ~ const (10) 

which explicitly shows that D((p q ) is inversely proportional to the length of 
the cascade. 

We have performed numerical simulations of the a-model in order to ob- 
tain the feeling to what extent these theoretical prejudices are realized in 
practice. In Figs. 1, 2 the histograms of the values of intermittency exponent 
<£> 2 , y?3 are plotted for 5000 generated cascades with 6 and 10 steps. These 
numbers were determined for each event using the method described in [|TJ 
(and applied there to the JACEE event ||). One sees that both the average 
value and the dispersion depend on number of steps in the cascade. For 
small number of steps the average value obtained from simulation is smaller 
than the "true" value given by Eq. (7). For 10 steps, however, the simulation 
gives the average rather close to the theoretical result. The dispersion of 
the distribution estimated directly from the observed peak, decreases with 
the number of cascade steps following well the 1/n rule of the Eq.(10). Its 
numerical value as a function of the cascade length is presented in Tables 1, 
2 for 2 different sets of cascade parameters a, b. The dispersion is relatively 
small, and it allows to distinguish between the cascades with different pa- 
rameters (Figs.l, 2). 
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Table 1. Intermittency exponents and their dispersions for a = 0.8, b — 1.1 
and n — 5, ... ,10 cascade steps 





theor. 


5 


6 


7 


8 


9 


10 


V2 = I0" 2 x 


2.9 


2.4 ±0.9 


2.5 ±0.8 


2.6 ±0.7 


2.7 ±0.6 


2.7 ±0.6 


2.7 ±0.5 


^3 = 10" 2 X 


8.2 


6.9 ±2.6 


7.1 ±2.3 


6.6 ±2.0 


7.7 ± 1.7 


7.7 ±1.6 


7.8 ±1.5 



Table 2. Intermittency exponents and their dispersions for a = 0.5, b = 1.5 
and n — 5, ... ,10 cascade steps 





theor. 


5 


6 


7 


8 


9 


10 




3.2 


2.4 ± 1.0 


2.5 ±1.0 


2.5 ±0.8 


2.7 ±0.7 


2.7 ±0.6 


2.8 ±0.6 


<p 3 = lO^x 


8.1 


5.9 ±2.3 


6.1 ±2.2 


6.4 ±2.1 


6.7 ± 1.8 


6.7 ± 1.6 


6.8 ±1.6 



4. Our conclusions can be summarized as follows : 

(a) the average value of the intermittency exponent obtained from our 
analysis is fairly close to the "theoretical" value. 

(b) the dispersion of the distribution is inversely proportional to the 
length of the cascade. It is found to be relatively small. This allows to 
distinguish between cascades with reasonably different parameters. 
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Figure captions 



Fig. 1 Distribution of the intermittency exponent ip 2 as determined in indi- 
vidual events generated from the a— model. 5000 events with a = 0.8, b — 1.1 
(a) and 5000 events with a = 0.5, b = 1.5 (b) were used. Cases (a) and (b) 
are plotted in two different scales. Histogram for the case (a) is multiplied 
by 10 _1 . 

Solid line and dots : 10 cascade steps, dashed line and crosses : 6 cascade 
steps. 

Fig. 2 Distribution of the intermittency exponent </? 3 . Other details as in 
Fig. 1. 
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